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9- Write a MATLAB code to design an ANFIS model to approximate the following 
function: 


fix ) = -lx - x 2 , where x e [-10 , 10]. 

Show the best values of the optimized parameters after leaning or training is done. 
To build the ANFIS model, consider the following Specs: 

> Use 60 training data pairs of input and output. 

> Use 3 MFs of triangular type for the input x. 


> Use the BP algorithm as a learning (or optimization) algorithm. 



> The no. of iterations (or epochs) in the training process = 100 



> The total No. of rules = 3 rule 

> The triangular MF has 3 control parameters. 

> No. of premise parameters = No. of control parameters for MF x total No. of MFs 

= 3x3 = 9 parameter 

> No. of consequent parameters = (No. of inputs + 1) x No. of rules 

= (1 + 1) x 3 = 6 parameter 

> The total No. of ANFIS parameters 

= No. of premise parameters + No. of consequent parameters 
= 9 + 6=15 parameter 

% First, we randomly obtain 60 training data pairs of input and output 
% use the following form of random distribution 

% Xi = Xi L + ( Xi U - Xi L ).*rand (No. of required points , 1); 

» xl = -10 + (10+10). *rand(60,l); % xl is a column vector with dimension 60x1 
» y = -2.*xl- xl.*xl; % y is a column vector with dimension 60x1 
» x = [xl y]; % the training data set 

» n = 60; % no. of samples of the training data set 

» numEpochs = 100; 
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» numMFs = 3; 

» mfType = ' trimf % we can use 'gaussmf ' or ' trapmf ' or 'gbellmf ' 

% build ANFIS model with the default initialization of parameters 

% Save in itial ANFIS model as fismatl 

» fismatl= genfisl(x, numMFs, mfType); 

% Or we can use: fismatl=genfisl(x,[3 5],char('trimf ','gaussmf ')); 

% Train the initial ANFIS model over the training data set x 

% Save the trained ANFIS model as fismat2 

» [fismat2,trn_mse,tst_mse] = anfis(x, fismatl, numEpochs,NaN,x,0); 

% NaN represent the default values of training options like training epoch number, 

% training error goal, initial step size, step size decrease rate and step size increase rate 

% trn_mse is avector contains the values of mean square errors during the epochs of 
training process 

% tst_mse is avector contains the values of mean square errors during testing process 
% 0 is used for the BP learning algorithm and 1 for the HL learning algorithm 
» tm_out = evalfis(x(:,l),fismat2); % the values of outputs for the trained ANFIS 
» [x(:,2) trn_out] % view the desired output and the trained ANFIS output 
» AMSE = mean(trn_mse); % the average of mean square errors over the 100 epochs. 
» epoch = ImumEpochs; 

» plot(epoch, trn_mse) % curve for mean square error vs. epochs 
» yy = [1 : n ]; % Data Set Index 

» plot(yy,x(:,2),Y,yy,trn_out,'o') 

10- In problem 2, state which of the following steps improve the approximation 
results for ANFIS: 

> Change the No. of MFs for the input x to 5. 

> Change the type of MFs to Gaussian. 

> Chang the learning algorithm to HL. 

> The total No. of rules = 5 rule 

> The Gaussian MF has 2 control parameters. 

> No. of premise parameters = No. of control parameters for MF x total No. of MFs 
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= 2x5 = 10 parameter 

> No. of consequent parameters = (No. of inputs + 1) x No. of rules 

= (1 + 1) x 5 = 10 parameter 

> The total No. of ANFIS parameters 

= No. of premise parameters + No. of consequent parameters 
= 10+10 = 20 parameter 

» xl= - 1 0+( 10+10). *rand(60, 1 ) ; 

» y= -2.*xl-xl.*xl; 

» x=[xl y]; 

» n=60; 

» numEpochs=100; 

» numMFs=5; 

» mfType='gaussmf ; 

» fismat 1 =genfis 1 (x,numMFs ,mfT ype) ; 

» [fismat2,trn_mse,tst_mse]= anfis(x,fismatl,numEpochs,NaN,x,l); 

» tm_out= evalfis(x(:,l),fismat2); 

» [x(:,2) trn_out] 

» AMSE= mean(trn_mse); 

» epoch= EnumEpochs; 

» plot(epoch,trn_mse) 

» yy=[l:n]; 

» plot(yy,x(:,2),'x',yy,trn_out,'o') 

11. Simulate the behavior of the system (the system output is Jt): 

x=- 0.6 x +f(x) 

Where / (jf) = -2x - X z is the same as in problem 8. 

First, build your simulink model with the exact expression for f(x ) then, replace f (x) 
by the ANFIS fuzzy model that you reached to in problem 8 and compare the 
results. To run the simulink model use the following initial values: 

x(0) = - 0.2 and x(0)= 0.5 

REPORT 


Fuzzy Logic Control (FLC) 


